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Abstract 

We suggest a generalization of vector calculus for the case of non-integer dimensional 
space. The first and second orders operations such as gradient, divergence, the scalar and 
vector Laplace operators for non-integer dimensional space are defined. For simplification 
we consider scalar and vector fields that are independent of angles. We formulate a 
generalization of vector calculus for rotationally covariant scalar and vector functions. 

This generalization allows us to describe fractal media and materials in the framework 
of continuum models with non-integer dimensional space. As examples of application 
of the suggested calculus, we consider elasticity of fractal materials (fractal hollow ball 
and fractal cylindrical pipe with pressure inside and outside), steady distribution of heat 
in fractal media, electric field of fractal charged cylinder. We solve the correspondent 
equations for non-integer dimensional space models. 


1 Introduction 

In general we can assume that space and space-time dimensions are D , which need not be an 
integer. Non-integer dimensional spaces and method of dimensional continuation are initially 
emerged in statistical mechanics and quantum held theory. Non-integer dimension D = 4 — 
e of space-time and e-expansion are actively used in the theory of critical phenomena and 
phase transitions in statistical physics (for example, see [I], 2j). Integration over non-integer 
dimensional spaces is used in the dimensional regularization method as a powerful tool to obtain 
exact results without ultraviolet divergences in quantum held theory ME]. In quantum 
theory, the divergences are parameterized as quantities with coefficients e -1 = (4 — D )~ l , and 
then these divergences can be removed by renormalization to obtain finite physical values. 

The axioms for integrals in non-integer dimensional space are suggested by Wilson in [6J. 
These properties are natural and necessary in applications in different areas |5j. Theory of 
integration in non-integer dimensional spaces has been suggested in jTj [5, ESj- Stillinger intro¬ 
duces [7\ a mathematical basis of integration on spaces with non-integer dimensions. In [7] 
has been suggested a generalization of the Laplace operator for non-integer dimensional spaces 
also. In the book by Collins [S] the integration in non-integer dimensional spaces is formulated 
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for rotationally covariant functions. The product measure method, which is suggested in [9], 
and the Stillinger’s approach [7] are extended by Palmer and Stavrinou [8] to multiple variables 
and different degrees of confinement in orthogonal directions. In the paper [8] extensions of 
integration and scalar Laplace operator for non-integer dimensional spaces are suggested. 

The scalar Laplace operators, which are suggested in mi for non-integer dimensional 
spaces, have a wide application in physics and mechanics. Non-integer dimensional space has 
successfully been used as an effective physical description. The Stillinger’s form of Laplacian 
first applied by He [TO, TTJ [T2J T3j, where the Schrodinger equation in non-integer dimensional 
space is used and the real confining structure is replaced by an effective space, such that the 
measure of the anisotropy or confinement is given by the non-integer dimension. Non-integer 
dimensions is used by Thilagam to describe stark shifts of excitonic complexes in quantum 
wells pH], exciton-phonon interaction in fractional dimensional space [15], and blocking effects 
in quantum wells [16 j. The non-integer dimensional space approach is used by Matos-Abiague 
mmmmmmm- to describe momentum operators for quantum systems and Bose-like 
oscillator in non-integer dimensional space, the polaron effect in quantum wells. Quantum me¬ 
chanical models with non-integer (fractional) dimensional space has been suggested by Palmer 
and Stavrinou |8], Lohe and Thilagam [23] • The non-integer dimensional space approach is used 
to describe algebraic properties of Weyl-ordered polynomials for the momentum and position 
operators [251 [26] and the correspondent coherent states [27]. The Stillinger’s form of Laplacian 
has been applied to the Schrodinger equation in non-integer dimensional space by Eid, Mus- 
lih, Baleanu, Rabei in [281122]; Muslih and Agrawal PH 151] . by Calcagni, Nardelli, Scalisi in 
[32]. The fractional Schrodinger equation with non-integer dimensions is considered by Martins, 
Ribeiro, Evangelista, Silva, Lenzi in [33] and by Sandev, Petreska, Lenzi [33]. Recent progress 
in non-integer dimensional space approach includes the description of the scalar field on non¬ 
integer dimensional spaces by Trinchero [35], the fractional diffusion equation in non-integer 
dimensional space and its solutions are suggested in [36]. The gravity in fractional dimensional 
space is described by Sadallah, Muslih, Baleanu in [37, [35], and by Calcagni in [32, HE SI] - 
The electromagnetic fields in non-integer dimensional space are considered in [ 321135 , 33] and 
[35] S61 SZl HU 111]. 

Unfortunately, the basic articles Ei proposed only the second order differential operators 
for scalar fields in the form of the scalar Laplacian in the non-integer dimensional space. The 
first order operators such as gradient, divergence, curl operators, and the vector Laplacian are 
not considered in Ei- In the book [39] (see also (35[ SSI HO SB]), the gradient, divergence, 
and curl operators are suggested only as approximations of the square of the Palmer-Stavrinou 
form of Laplace operator. Consideration only the scalar Laplacian in non-integer dimensional 
space approach greatly restricts us in application of continuum models with non-integer dimen¬ 
sional space for fractal media and material. For example, we cannot use the Stillinger’s form 
of Laplacian for displacement vector field u(r, t) in elasticity and thermoelasticity theories. We 
cannot consider equations for the electric field E(r, t) and the magnetic fields B(r, £) for elec¬ 
tromagnetic theory of fractal media by using continuum models with non-integer dimensional 
space. 
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In this paper, we propose a vector calculus for non-integer dimensional space and we define 
the first and second orders differential vector operations such as gradient, divergence, the scalar 
and vector Laplace operators for non-integer dimensional space. For simplification we consider 
rotationally covariant scalar and vector functions that are independent of angles. In order to 
derive the vector differential operators in non-integer dimensional space we use the method of 
analytic continuation in dimension. For this aim we get equations for these differential opera¬ 
tors for rotationally covariant functions in R n for arbitrary integer n to highlight the explicit 
relations with dimension n. Then the vector differential operators for non-integer dimension D 
are defined by analytic continuation from integer dimensions n to non-integer D. These pro¬ 
posed equations allows us to reduce D -dimensional vector differentiations to usual derivatives 
with respect to one variable r = |r|. It allows us to reduce differential equations in non-integer 
dimensional space to ordinary differential equations with respect to r. The proposed opera¬ 
tors allows us to describe fractal materials and media in the framework of continuum models 
with non-integer dimensional spaces. In order to give examples of the possible applications, 
we consider continuum models of fractal media and materials in the elasticity theory in the 
heat theory, and in the theory of electric fields. The correspondent equations for non-integer 
dimensional space are solved. 


2 Fractal media 

The cornerstone of fractal media is the non-integer dimension [5D] such as mass or charge 
dimensions mm- In general, fractal media and materials can be treated with three different 
approaches: (1) Using the methods of’’Analysis on fractals” [53115T1155115611571158] it is possible 
to describe fractal materials; (2) To describe fractal media we can apply fractional-integral 
continuous models suggested in [59l [601 EH EH E2] (see also [63l [611 ESI [66l [671 EH El] and [701 
m m m m ca ebd. In this case we use integrations of non-integer orders and two different 
notions such as density of states and distribution function [S2j; (3) Fractal materials can be 
described by using the theory of integration and differentiation for a non-integer dimensional 
space 

The first approach, which is based on the use of analysis on fractal sets, is the most stringent 
possible method to describe idealized fractal media. Unfortunately, it has two lacks. Firstly, 
a possibility of application of the analysis on fractals to solve differential equations for real 
problems of fractal material is very limited due to weak development of this area of mathematics 
at this moment. Secondly, fractal materials and media cannot be described as fractals. The 
main property of the fractal is non-integer Hausdorff dimension m that should be observed on 
all scales. The fractal structure of real media cannot be observed on all scales from the infinitely 
small to the infinitely large sizes. Materials may have a fractal structure only for scales from 
the characteristic size of atoms or molecules of fractal media up to size of investigated sample 
of material. 

The second approach, which is based on the use of fractional integration in integer dimen- 
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sional spaces, can give adequate models to describe fractal media. The main disadvantage of 
the fractional-integral continuum models is the existence of various types of fractional integrals, 
which led to the arbitrariness in the choice of the correspondent densities of states. 

In this paper, we consider the third approach used the non-integer dimensional spaces. One 
of the advantages of this approach is a possibility to avoid the arbitrary choice of densities 
of states. In addition, we also suggest a generalization of the vector calculus to the case of 
non-integer dimensional space. It allows us to use different continuum models of fractal media 
and materials in the framework of non-integer dimensional space approach. 

Real fractal materials can be characterized by the relation between the mass Mo(W) of a ball 
region W of fractal medium, and the radius R of this ball in the form M D (W) = M 0 (R I / Rq ) d , 
R/Rq 3> 1, where Rq is the characteristic size of fractal medium such as a minimal scale of 
self-similarity of a considered fractal medium. The number D is called the mass dimension. 
The parameter D, does not depend on the shape of the region W, or on whether the packing 
of sphere of radius Rq is close packing, a random packing or a porous packing with a uniform 
distribution of holes. As a result, fractal materials can be considered as media with non¬ 
integer mass dimensions. Although, the non-integer dimension does not reflect completely the 
geometric and dynamic properties of the fractal media, it nevertheless allows us to get important 
conclusions about the behavior of these media. As it will be shown in the next section, the 
power law Mjj(W) ~ R D can be naturally derived by using the integration over non-integer 
dimensional space, where the space dimension is equal to the mass dimension of fractal media. 

In order to describe fractal media by continuum models with non-integer dimensional space, 
we should use the concepts of density of states c n (D, r) and distribution function p(r). The 
density of states describes how closely packed permitted places (states) in the space M”, where 
the fractal medium is distributed. The expression c n (D,r)dV n is equal to the number of per¬ 
mitted places (states) between V n and V n + dV n in M n . The distribution function describes a 
distribution of physical values such as mass, electric charge, number of particles on a set of 
permitted places (possible states). In general, the concepts of density of states and distribution 
function are different. We cannot reduce all properties of the fractal media to the distribution 
function only, and we should use the concepts of density of states to characterize how closely 
packet permitted states places in the media. 

The most important property of fractal medium is the fractality, which means that the mass 
M d (W) of this medium in any region W C M n increases more slowly than the n-dimensional 
volume V n (W) of this region. For the ball region of the fractal medium, this property can be 
described by the power law Md(W) ~ R D , where R is the radius of the ball, and D is the mass 
dimension. 

Another important property of some fractal media is homogeneity. Fractal medium is called 
homogeneous if the power law Md(W) ~ R D does not depend on the translation of the region. 
The homogeneity property of the fractal medium means that two regions W\ and IT 2 with the 
equal volumes V n {W\) = 14 (RA) have equal masses Ad^iWi) = Md(W 2 ). 

To adequately describe the fractal media by continuum models with non-integer dimensional 
spaces, the following two requirements must be satisfied. 
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• In the continuum models the mass density of homogeneous fractal medium should be 
described by the constant distribution function p(r) = po = const. Then equations with 
constant density should describe the homogeneous media, i.e., the conditions p( r) = const 
and V n {Wi) = V n {W 2 ) should lead to the relation Md{W\) = Md{W 2 ). 

• In the continuum models the mass of the ball region W of fractal homogeneous medium 
should be described by a power law relation M ~ R D , where 0 < D < 3, and R is the 
radius of the ball. Then the conditions V n {W\) = \ n V n {W 2 ) and p( r) = const, should 
lead to the relation Mu{W \) = \ D Md(W 2 ). 

These requirements cannot be realized if the mass of fractal medium is described by inte¬ 
gration of integer order over the integer dimensional space without using the concept of density 
of states c n (D, r). In order to realize these requirements we propose to use the integration and 
differentiation in non-integer dimensional spaces. In this case we can use the equation 

dM D (W) = p( r) dV D ( r, n), (1) 

where p( r) is a distribution function, and the density of states c n {D, r) in M n is chosen such 
that 

dV D ( r, n) =Cn{D , r) dV n 

describes the number of permitted states in d,V n . For different values of n G {1,2,3} we can 
use the notations 


dV D = c 3 (D,t) dV 3 , dS d = c 2 (d,r) dS 2 , dip = ci(/3,r) dh (2) 

to describe fractal media in the spaces R", where these media are distributed. For simplification, 
we also use the notation d D r instead of dVo(r, n). The form of function c n (D, r) is defined by 
the properties of considered fractal medium. The symmetry of the density of states c n (D, r) 
is dictated by the symmetry properties of the described fractal medium, but in any cases it 
should be a function of power-law type to adequately reflect a scaling property (fractality) of 
described fractal medium. To simplify the analysis in this article we will consider only isotropic 
fractal media with densities of states that are independent of angles. 

In the continuum models of fractal media, it is convenient to work in the dimensionless space 
variables x/R 0 — > x, y/Ro — > x, z/R 0 — > x, r /R 0 —» r, that yields dimensionless integration 
and dimensionless differentiation in non-integer dimensional space. In this case the physical 
and mechanical quantities of fractal media have correct physical dimensions. 

3 Integration over non-integer dimensional space 

The integral for all non-integer values of D is defined by continuation in D Si. Let us 
give properties must we impose on a functional of /(r) in order to regard it as D -dimensional 
integration. The following properties or axioms [6] for integrals in D -dimensional space are 
natural and necessary in applications [5]: 


5 


1 . 


Linearity: 


J («/i(r) + &/ 2 (r)) d D r 


= a 



(3) 


where a and b are arbitrary real numbers, and d D r = dVo(r, n) represents the volume 
element in the non-integer dimensional space.. 


2. Translational invariance: 


J f(r + r 0 ) d D r 



(4) 


for any vector ro- 
3. Scaling property: 



—D 


J f(r)d D r 


(5) 


for any positive A. 

Linearity is true of any integration, while translation and rotation invariance are basic 
properties of a Euclidean space. The scaling property embodies the D- dimensionality. Not 
only are the above three axioms necessary, but they also ensure that integration is unique, 
aside from an overall normalization [6], 

These properties must be imposed on a functional of /(r) in order to regard it as D- 
dimensional integrations [5]. These properties are natural and are necessary in application of 
dimensional regularization to quantum field theory (see section 4 in [5]). 

A function /(r) that we integrate could in principle be any function of the components of 
its vector argument r. However, we do not, a priori, know the meaning of the components of, 
say, a vector in non-integer dimensions. In this paper, we will work with rotationally covariant 
functions. So we will assume that / is a scalar or vector function only of scalar products of 
vectors or of length of vectors. For example, in the elasticity theory, we consider the case, 
where the displacement vector u(r), is independent of the angles u(r) = u(r), where r = |r|. 

The integral defined in equation (J6]) satisfies the properties ([3]) - ([5]) . 

The .D-dimensional integration for scalar functions /( r) = /(|r|) can be defined in terms of 
ordinary integration by the expression 




where we can use 



(7) 


For integer D = n, equation (J7J) gives the well-known area S n _i of (n — l)-sphere with unit 
radius. 
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As a result, we have [5] the explicit definition of the continuation of integration from integer 
n to arbitrary fractional D in the form 

r r°° 

j d D r /(|r|) = —-- J dr r D ~ 1 f(r). (8) 

This equation reduced D-dimensional integration to ordinary integration. Therefore the lin¬ 
earity and translation invariance follow from linearity and translation invariance of ordinary 
integration. The scaling and rotation covariance are explicit properties of the definition. 

Let us give some examples of an application of equation (H). For the function 


/(r 2 ) = 


r 2 + a 


r 2 + 6’ 

where a and b are real numbers. The integral for (J9]) can be explicitly computed 


D r + a 


r v+» 

The other example is the integral 


= {irb) D/2 (a/6 - 1) T(1 - D/2). 


d D r 


„2 a 


r ““ _ T(a + D/2) T(/3 a D/2) n /2 D+ 2 a- 2/3 

(r 2 + a 2 )? ~ T(D/2)T(/3) 


(9) 


( 10 ) 


( 11 ) 


where r = |r|. 

The mass of material in W is described by the integral 

M d (W)= [ p{v)d D r, 
Jw 


( 12 ) 


where r is dimensionless vector variable. For a ball with radius R and the density p(r) = po = 
const, we get the mass is defined by 


D /2 

M d (W ) =p 0 V D = 7 , P \ R d . 
> FO r(D/2 + l) 


(13) 


This equation defines the mass of the fractal homogeneous ball with volume V/> For D = 3, 
equation (TT3|) gives the well-known equation for mass of non-fractal ball M 3 = ( 4:p 0 7i/3)R 3 
because T(3/2) = y/n/2 and T(z + 1) = zT(z). 


4 Divergence for non-integer dimensional space 

Let us consider hollow ball Bo(Ri] R 2 ) with internal radius R\ and external radius R2 in 
non-integer dimensional space. The boundary dB D (Rp R 2 ) of this ball consists of two (D — 1)- 
dimensional spheres Sd-i(Ri) and S'z ?—1 (-R 2 )- 
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We assume that the vector held u(r) is radially directed and u(r) is not dependent on the 
angles, i.e., 

u(r) = u r (r) e r , (14) 

where e r = r/|r| and r = |r|. We can dehned a hux of the vector held u(r) across a surface 
OBd^Ri] R 2 ) by using the integration in non-integer dimensional space suggested in [5]. A hux 
of the vector held u(r) across a (D — l)-dimensional surface So -1 is the integral 

®u(Sd-i) = [ (u,dS D _i), (15) 

Js D _! 


where dS D _i = e r dSA-i. 

The volume of B D {R ,\; R 2 ) is equal to 


V(B d (R 1 -R 2 )) = 


7r 


D/2 


(i? 2 D - R?' 


T{D/ 2 + 1) 


(16) 


An exact expression for dependence of the divergence operator on the non-integer dimension 
D and the vector held u(r) = u r (r ) e r can be derived by the equation 


Divf u = 


lim 


^ u {dB D {R i; R2)) 


V(B d (i?i;i? 2 ))—>-0 V(B D (R l ]R 2 )) 
Here the hux of the vector held can be represented by 


(17) 


l>„(9 B d (RuR 2 ))= / (u, dS D _ 1 ) 

J <Sd_i(-R 2 ) 


' Sd-i(Ri) 


(u,dS D _i). 


(18) 


Using ([HD , we get 


ldB D (R v ,R 2 ) 


(u, dS D -i) = 


[ (u, dS D -i) — [ (u, 

Js D _ 1 (R 2 ) J Sd-i(-Ri) 


dSo.,) = 


= u r (R 2 ) / dS D _i- u^R^ / dS D _i 

Js D _!(i? 2 ) JSd-^R!) 


2 7r D / 2 


(u r (R 2 ) R/- 1 - u r (Ri) R/-‘y 


TiD/2) 

where 5 ' D -i = (e r ,S D _i). 

For D = 3, we have T(3/2) = (1/2) T( 1/2) = y/n/2, and 


(19) 


ldB D (Rv,R 2 ) 


(u, dSo-i ) = 47t (u r (R 2 ) i ?2 — u r (R\) R\ 


( 20 ) 
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For an infinitely thin hollow ball (thin spherical shell) with R\ — r and R 2 = r + Ar, we 
can use 

,d . 


u $ u (dB D (r-,r + Ar)) 

Div r u = Inn Tr . . - - 

Ar->o V{B£,{r\r + Ar)) 


to derive an expression for the divergence for our case. 
The volume of the ball Bn{r\ r + Ar) is 

D/2 


7T"' 1 ' / \ 

1/(5 D (r; r + Ar)) = p — yy (V + ~~ r D ) = 


71 


D/2 


/ \ 97 B / \ 

r D ~ l Ar + 0((Ar) 2 )J = ——— (rAr + 0((Ar) 2 ) J. (22) 


r(D /2 + l)V vv r(D/ 2 ) 

The flux of the vector field u = u r {r ) e r across a surface dB D (Ri, R 2 ) is given by 


( 21 ) 


2ir D ' 2 

T(B/2) 


4> u (SB D (r|r +Ar)) = / (u,<iS D _ 1 ) - / (u,<iS D _ 1 ) = 

J S D _i(r+Ar) 7s_D_i(r) 


Ur(r)- 


2tt d / 2 

"@j 

<9w r (r) 


^u r (r + Ar) (r + Ar) D 1 «- u r {r)r D 


<9r 


Ar+0((Ar) 2 ) r D 1 + (£) — 1) r 15 2 Ar+0((Ar) 2 ) — u r (r) r D ^ 


pT /2) ( rD 1 ^5^ + ( B - 1) rD 2 “r(r) + C(Ar)) Ar. 


r(£>/ 2 ) 

Using (1221) and (l23lh we get 

<& u (dB D (r- : r + Ar)) du r (r) D — 1 


V(B D (r] r + Ar)) <9r 


+ 


w r (r) + O(Ar). 


(23) 


(24) 


As a result, the divergence for the vector field u = u r (r) e r in non-integer dimensional space 
has the form 


Divf* u = lim , 

Ar— >0 V(Bd{t] r + Ar)) 


§ u (dB D (r; r + Ar)) du r (r ) D — 1 


+ 


■u r (r). 


dr r 

The Gauss theorem in non-integer dimensional space can be written in the form 


(25) 


/ Divf u d D r = / (u, dS^.i), 

'B D (RyR 2 ) JdB D {RyR 2 ) 


(26) 


where we assume that the dimension .D of the region of fractal materials and the dimension d 
of boundary of this region are related by the equation d = D — 1. 
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5 Vector differential operators in non-integer dimensional 
space 


We would like to derive equations for vector differential operators in non-integer dimensional 
space. For this aim we should have equations for these differential operators for rotationally 
covariant functions in the spherical coordinates in M” for arbitrary n to highlight the explicit 
relations with dimension n. Then the vector differential operators for non-integer dimension D 
can be defined by analytic continuation in dimension from integer n to non-integer D. 

To simplify we will consider only scalar fields p and vector fields u that are independent of 
angles 

<p(r) = <p(r), u(r) = u(r) = u r e r , 

where e r = r/r, r = |r| and u r = u r (r) is the radial component of u. We will work with 
rotationally covariant functions only. This simplification is analogous to the simplification for 
definition of integration over non-integer dimensional space suggested in [5]. 

Let us give equations for differential operators for functions u = u r (r) e r and p = p(r) in 
the spherical coordinates in M n for arbitrary n. 

The divergence in integer dimensional space M n for the vector held u = u(r) e r is 


,. du r n — 1 

chv u = —-1- u r . 

or r 

The gradient in integer dimensional space M rt for the scalar held p = <p(r) is 


(27) 


grad ip = rE- e r 
or 


(28) 


The scalar Laplacian in integer dimensional space W 1 for the scalar held p = p(r) is 


A ip(r) = divgrad^(r) 


d 2 p n — 1 dp 
dr 2 r dr 


(29) 


The vector Laplacian [78] in integer dimensional space M rt for the vector held u = u(r) e r is 


Au(r) = graddivu(r) 


rd 2 u r 
V dr 2 ~*~ 


7i — 1 du r 
r dr 



(30) 


As a result, we have equations of differential operators in K n for continuation from integer 
n to arbitrary non-integer D in the following forms. 

The divergence in non-integer dimensional space for the vector held u = u(r) is 


Divf u 



D - 1 
- u r . 

r 


(31) 


The gradient in non-integer dimensional space for the scalar held p = p{r) is 


n i D di f 

Grad r p = — e r . 
dr 


(32) 


10 












( 33 ) 


The curl operator for the vector field u = u(r) is equal to zero 

Curl” u = 0. 


The scalar Laplacian in non-integer dimensional space for the scalar field <p = ip(r) is 


A = Div" Gradf ip = 


d 2 (p D 


+ 


1 dip 
dr 


2 ry* 

The vector Laplacian in non-integer dimensional space for the vector field u = u{r) e r is 

1 du r D — 1 


(34) 


v Af u — Gradr Div^ u — 


D 


d 2 u r D 


2 ry* Qry> ry* ~ 


Un 


G r . 


(35) 


Let us consider a case of axial symmetry of the fractal material, where the fields <p(r) and 
u(r) = u r (r) e r are also axially symmetric. Let Z-axis be directed along the axis of symmetry. 
Therefore it is convenient to use a cylindrical coordinate system. Equations for differential 
vector operations for cylindrical symmetry case have the following forms. 

The divergence in non-integer dimensional space for the vector field u = u(r) is 

2 (36) 


n du r D 
Div” u = 4- 


■ U r - 


dr r 

The gradient in non-integer dimensional space for the scalar field ip = <p(r) is 


Grad” <p = Af- e r . 

dr 

The scalar Laplacian in non-integer dimensional space for the scalar field ip = <p(r) is 


dip 


(37) 


c.n T^- ^D D — 2dip 

(p = Divy Grad” <p = ^^7 + 


dr 2 r dr 

The vector Laplacian in non-integer dimensional space for the vector field u = u(r) e r is 


(38) 


v a d 


Af u = Grad” Div, 


,D 


U = ( 


d 2 u r D 

,- - 4- 

V dr 2 


2 du r 


D — 2 


dr 


u, 


(39) 


For D — 3 equations ( 1311 ) - ( 1391 ) give the well-known expressions for the gradient, divergence, 
scalar Laplacian and vector Laplacian in M 3 for fields (p = < p(r ) and u(r) = u r (r) e r . 

It is easy to generalize these equations for the case <p = <p(r, z ) and u(r, z) = u r (r, z ) e r + 
u r (r, z) e z . In this case the curl operator for u(r, z) is different from zero, and 


Curl”., u = 


du r du , 


dz dr 


(40) 


The suggested operators allow us to reduce D- dimensional vector differentiations (1TB - (ITT]) 
and (ITT]) - flTTj) to derivatives with respect to r = |r|. It allows us to reduce partial differential 
equations for fields in non-integer dimensional space to ordinary differential equations with 
respect to r. 
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6 Stillinger’s Laplacian for non-integer dimensional space 


In the paper [7], the integration in a non-integer dimensional space is described by using the 
equation 


27r( D_1 )/ 2 f°° 

d = r((D- l)/2) / 

where Gpr = dVoiv^n) represents the volume element in the non-integer dimensional space. 
Using (14B with ip(r,8 ) = 1, and 



drr 


D -1 


dd <p(r,6) sin D 2 0 , 


(41) 


sin® 2 6 = 


7T 1 / 2 r(D/2- 1) 

rp/ 2 ) 


(42) 


we get 

vr^/2 

which is the volume of a D -dimensional ball with radius R. 

Using the analytic continuation of Gaussian integrals the scalar Laplace operator for non¬ 
integer dimensional space has been suggested in (7J. For a function ip = <p(r, 6) of radial distance 
r and related angle 6 measured relative to an axis passing through the origin, the Laplacian in 
non-integer dimensional space proposed by Stillinger [7] is 



St^D 


1 d ( D _ 1 d\ 1 d ( . D _ 2 d \ 

r^dr V fr) + r 2 sin D - 2 0d# \ 9 dd) ’ 


(44) 


where D is the dimension of space (0 < D < 3), and the variables r > 0, 0 < 9 < n. Note that 
(^A D) 2 ^ st A 2D . If the function depends on radial distance r only ((y9 = then 


St A D (p(r) 


1 9 ( D _ x d(p(r) \ 

r D ~ i dr \ dr J 


d 2 <p(r) D — 1 d(p(r) 

dr 2 r dr 


(45) 


It is easy to see that the Stillinger’s form of Laplacian st A D for radial scalar functions </?(r) = 
<p(r) coincides with the suggested scalar Laplacian s A® for this function, 

st A D (p(r ) = S A D tp(r). (46) 


The Stillinger’s Laplacian can be applied only for scalar fields and it cannot be used to describe 
vector fields u = u r (r) e r because Stillinger’s Laplacian for D = 3 is not equal to the usual 
vector Laplacian for R 3 , 


St A 3 u (r) ± Au(r) 


/ d 2 u r 
V dr 2 


+ 


2 du r 


r dr 



6 r . 


(47) 


For the vector holds u = u r (r ) e r , we should use the vector Laplace operators (15^1) and 

Note that the gradient, divergence, curl operator and vector Laplacian are not considered 
in Stillinger’s paper |7J 
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7 Non-integer dimensional space for complex fractal me¬ 
dia with d D — 1 

In general, the dimension D of the ball region Bo of fractal materials and the dimension d of 
boundary Sd = dB D of this region are not related by the equation d — D — 1, i.e., 

dim(dB D ) 7 ^ din a(B D ) — 1, (48) 

where dim (Bd) = D. We assume that dimension of the boundary Sd = dBo is 

dim(5'd) = d. (49) 

Considering an infinitely thin hollow ball Bo with R\ — r and R% = r + Ar, we can use 


Div?’ d u = lim 


®u{Sd{r-,r + A r)) 


Ar-K) V ( Bo{r ; r + Ar)) 

in order to derive an expression for the divergence for the case d 7 ^ D — 1. 
The flux of the vector held u = uJr) e r across a surface Sd is 


(50) 


®u(S<i(r;r + Ar)) = 


2 x ( d+1 )/ 2 

r«d+1)/2) 


Using 


2 7r ( rf + 1 )/ 2 
= r((d + l)/ 2 ) 

and (I5TT) . we get 


du r (r ) 
dr 


(u r (r + Ar) (r + A r) d — u r (r ) r d ^j = 
+ dr d ~ 1 u r {r) + O(Ar) j Ar. 


(51) 


u(*S'd(r; r + Ar)) ^~^T(D/2) , 1 du r (r) J_ x 

(5 D (r;r + Ar)) T((d + l)/2) Vr^" 1 ^ dr + j 


$ 
V(B_ 


Using the limit Ar —> 0, we get the divergence operator for fractal media with d D — 1 in 
the form 

n (d+i- d )/2 r( D /2) ( 1 du r [r ) d ( d\ 

( 53 ) 


Div?’ d u = 


T((d+l)/ 2 ) \r D ~ v 


)—1—d 


dr 


+ 


„D-d r 


u r (r. 


We can define the parameter 

a r = D — d, (54) 

that can be interpreted as a dimension of medium along the radial direction. Using (l54]h 
equation (1531) can be rewritten in the form 

Div?- 1 u = AA a -> /2 > ( -A AM + ± <r) ) . (55) 

T((d + l)/ 2 ) \r“ r_1 dr r Qr / 

This is divergence operator for non-integer dimensional continuum models of fractal materials 
with d 7 ^ D — 1. For ay = 1, i.e. d = D — 1, equations (l53|) and (l55|i give (jSIjl . 
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We can assume that the gradient for the scalar held <^(r) = p (r) may depend on the radial 
dimension a r in the form 


Gradf'V= r( ° r/ 2 ) 1 ^ (r) e,. 

/j^Oir / 2i iy*OLr X Qj. 

because expression (T56l) can be represented by the equation 

Grad?"V= lim ^ + , Ar) tk 

Ar— >0 V(B ar (r\ r + Ar)) 


where 


V ( B ar (r; r + Ar)) = —— ( r ar 1 Ar + 0((Ar) 2 ) 


(56) 


(57) 


(58) 


I>r/2) 

The presence of the factor of 2 in (157)) is due to the fact that for D — 1, r is integrated from 
—R to R , and when the limits are taken as 0 and R , one gets a factor of 2. 

For a r = 1, equation (1561) gives fj32|) . 

Using the operators (1561) and (1531) for the fields p = p(r) and u = u(r) e r , we can get the 
scalar and vector Laplace operators for the case d ^ D — 1 by the equation 


s A?- ■'V = Divf ■ 4 Gradf p, v A?' d u = GradfDivf u. 

Then the scalar Laplacian for d ^ D — 1 for the held p — p{r) is 

s D d = T((d + a r )/2) T(qy/2) / 1 d 2 p d + 1 - a r dp_ 

7 T ar ~ 1/2 T((d + 1)/2) \r 2ctr ~ 2 dr 2 r 2a r -i g r 


and the vector Laplacian for d ^ D — 1 for the held u = u{r) e r is 


V/\ D < d u — I a r)/^) r(oy/2) /_1 


7T ar_ 1/2 r((d + l)/2) \r 2 a r - 2 Q r 2 


+ d + 1 — ay <9u 


^2 ckt ’ 1 


<9r 


da r 

rv* OiOif 


u 7 


(59) 


(60) 


(61) 


The vector differential operators (l56jh (l53lh (1601) and (16TT) allow us to describe complex 
fractal materials with the boundary dimension of the regions d ^ D — 1. 


8 Applications in mechanics and physics 

In this section, we consider some applications of vector calculus in non-integer dimensional 
space to the elasticity theory, the heat processes, and electrodynamics. 

8.1 Elasticity of fractal material 

For homogenous and isotropic materials, the equation of linear elasticity [79], [80] for the dis¬ 
placement vector helds u = u(r, t ) has the form 

A grad div u + 2/rAu + f = p D 2 u, (62) 
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where A and /i are the Lame coefficients, and f is the external force density vector held. 

If the deformation in the material is described by u(r, t) = u(r, t ) e r , then equation (1621) has 
the form 

(A + 2/i) A u(r, t) + f (r,t) = p D 2 u(r, t ). (63) 

A generalization of equations (]63j) for fractal material in the framework of non-integer 
dimensional models has the form 

( A + 2/i) v Af u(r, t) + f (r, t) = p D 2 u(r , t ) , (64) 

where 1 A(r* is defined by (1351) . Equation (154)) describes dynamics of displacement vector for 
fractal materials. For static case, equation (164)) has the form 

V A? u(r) + (A + 2/i) -1 f(r) = 0, (65) 

where u = u r e r and f = f(r) e r . 

Let us consider some two problems for elasticity of fractal materials. 


Elasticity of fractal hollow ball with pressure inside and outside 


Let us determine the deformation of a hollow fractal ball with internal radius R\ and external 
radius i? 2 , with the pressure p\ inside and the pressure P 2 outside. 

We can use the spherical polar coordinates with the origin at the center of the ball. The 
displacement vector u is everywhere radial, and it is a function of r = |r| alone. Then the 
equilibrium equation for fractal ball is 


(A + 2/r) v A^ u(r) = 0, 

where u = u r e r . Using (135|) . we represent equation (1661) in the form 

d 2 u r (r) D — 1 du r (r ) D — 1 

H-^-o—ttr(r) = 0. 


dr 2 


dr 


The solution of (1671) is 


u{r ) = C\ r + U 2 r 1 D . 

The constants C\ and C 2 are determined from the boundary conditions for radial stress 


( 66 ) 

(67) 

( 68 ) 


(T rr (/?l) Pi, CTrr^Ro) P2- 


Using that the radial components of the stress is 


oy r (r) 


(2 p + A) 


du r 

dr 


+ A 


D - 1 
r 


we get 

,. -(piRS-PiR?) 

1 (2v + D\)(R°-R?y 


(69) 


(70) 

(71) 
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P2 ~ Pl 


(72) 


C 2 = 


2(l-D)n(R?-R?y 

Then the radial components of the stress is 


cr rr (r) = 


" (p2 7?” - P\ R\) , (P2 “Pi) (-R 1 -R 2 ) —D 


\D 


+ 


pD pO p D pfl 

it2 ./ti 1X9 il-| 


■ r 


(73) 


L 2 -^1 

The stress distribution in a ball with pressure p\ — p inside and p 2 = 0 outside is gives by 


CTrr IT = 


p R\ 


D 


pfl _ jdD 

k 2 


1 - 1 '* 




(74) 


The stress distribution in an infinite elastic medium with spherical cavity with radius R 
subjected to hydrostatic compression is 


Vrr(r) = -P I 1 - ( j 


D' 


(75) 


that can be obtained by putting R 1 = R, R 2 —» cxo, p l = 0 and p 2 = p in equation (175)1 . 


Elasticity of cylindrical fractal solid pipe with pressure inside and outside 

Let us consider the deformation of a fractal solid cylindrical pipe with internal radius Ri and 
external radius R 2 with pressure p\ a inside and pressure p 2 outside. We use the cylindrical 
coordinates with the z-axis along the axis of the pipe. When the pressure is uniform along 
the pipe, the deformation is a purely radial displacement u = u r {r) e r , where e r = r/r. The 
equation for the displacement u r (r) in fractal pine is 


d 2 u r (r) D — 2 du r 

2 fjp 


D -2 

—^— u, 


0 , 


where 0 < D < 3. If D = 3, we get the usual (non-fractal) case. 

The general solution of equation (1751) . where D ^ 1, D ^ 2, has the form 


(76) 


u r (r) = Ci r + C 2 r 2 ~ D . (77) 

Equations (1761) with D — 1 has the general solution 

u r {r) — Ci r + C 2 r ln(r). (78) 

For D = 2, equations (176|) has the solution 

u r (r) = Ci + C 2 r. (79) 


16 










Note that dimensions D = 1 or D = 2 of the fractal pipe material do not correspond to the 
distribution of matter along the line and surface. These dimensions describe a distribution of 
matter in 3-dimensional space (in the volume of pipe) such that the mass dimensions are equal 
to D. 

The constants C\ and C 2 are determined by boundary conditions. Using that pressure p\ a 
inside and pressure p 2 outside, we get the boundary condition in the form 


cr rr (i?x) pi, (J rr (R 2 ) p 2 . 


(80) 


Using (1771) and 


n_ o 

Orr = (2/1 + A)—- + A- u r = (2/i + A (D - 1)) C, - 2 p (D - 2 ) C 2 r l ~ D , (81) 

or r 


the boundary condition (15Uj) gives 

Ci = - 


Pi D - P2 R\ 


1 -D 


{2p + \{D-l)){R l 2 ~ D -R\~ D y 

r _ P2-P1 _ 

2 2 p(D — 2) (R\~ d — R\~ D )' 

The stress ior 2 < I) < 3 or 1 < D < 2 can represented in the form 


Pi R? 1 - P2 R 

(R ?- 1 - R?- 1 ) 




D—l 


P 2 ~ Pi 


(- R?~ 1 - 


Ri R 2 


D-l 


(82) 

(83) 

(84) 


For the boundary conditions o rr (R 2 ) = 0 and o rr (Ri) = —p, i.e. p 2 = 0 and pi = p for (1841) . 
we have the solution 




pRi 


D—l 


(fi?- 1 - fif- 1 ) 


1 - 1 '* 


D-r 


(85) 


This is the deformation of cylindrical pipe with a pressure p inside and no pressure outside. 
For D = 3, equation (1851) has the well-known form 


pR\ (, (R2 

{R\ - R\) ^ { r 

that describes the stress of non-fractal material of pipe. 


( 86 ) 


8.2 Heat equation for fractal materials 

The heat equation is 

diP i r ’ — ~ a Ap(r, t) = — q( r, t), 
ot c p p 


(87) 
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where <p(r,t) is the heat density of a medium, q(r,t) is the heat source density, and a is the 
thermal diffusivity a = k/c p p , where k is thermal conductivity, p is density, c p is specific heat 
capacity. 

A generalization of equation (l87|l for fractal material, has the form of the heat equation in 
the non-integer dimensional space 

d<P ^^ - a s A°ip(r,t) = — q(r,t), ( 88 ) 

dt c p p 


where we assume that the fields <p(r,t) and q(r,t ) are not depend on the angles. 
Using fHhfll . we get the following equations for the ball 


d 2 cp(r, t ) + D- 1 dq ?(r, t) + J_ , f) = 1 d( f ( r > t) 


dr' 2 r dr c v p 

For the pipe and cylinder, we get 

d 2 ip(r,t) D — 2dip(r,t ) 1 


+ 


dr 2 r dr c p p 

Steady states is described by the equation 


+- q(r,t) = 


a dt 


1 d(p(r,t ) 
a dt 


+ C - 1 9y(r) + J_ 
dr r dr c p p 

The general solution of equation (T9TT) is 

1 


ip(r) — Ci + C 2 r + 


c P p{D — 2) 


q{r)r u 1 dr — / g(r)rdr), (92) 


(89) 


(90) 


(91) 


where the constants C\ and Ci are determined by the boundary condition. For the constant 
heat source density q(r) = qo = const, equation (jS2]) has the form 


<p(r) = C\ + C 2 r 2 ~ D 


Qo 2 

-r 

2D c p p 


For D = 3, we get the well-known equation for non-fractal material. 


(93) 


8.3 Electric field of fractal charged infinite cylinder 

Let us consider a uniformly fractal charged infinite circular cylinder of radius R with a volume 
charge density p = const and non-integer dimension 2 < D < 3. Using the Poisson equation 
for scalar potential created by an infinite circular cylinder. We assume that the Z-axis is 
directed along the axis of the cylinder. Due to the axial symmetry of the charge distribution 
the potential is also axially symmetric. Therefore it is convenient to use a cylindrical coordinate 
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system. 

form 


The Poisson equation for scalar field cp(r) in non-integer dimensional space has the 


d 2 p D - 2 dp 
dr 2 r dr 


(0 <r<R), 
£o 


(94) 


ri‘ f D - 2 dtp _ 

dr 2 r dr 

The general solution of equations (1941) and (l95|l are 


(r > R). 


(95) 


r(r) = Ci + C 2 r 3 - D 


2 

pr 

2eo {D- 1) 


(0 < r < R), 


(96) 


if(r)=C 3 + C 4 r 3 - D (r > R), (97) 

where C\, C 2 , C 3 , C 4 are the integration constants, and 2 < D < 3. For the case D — 3, the 
general solution of equations (1941) and (j95|) has the well-known form 

2 

ip[r) = Ci + C 2 ln(r) - — (0 < r < R), (98) 

4e 0 

ip(r) = C 3 + C 4 ln(r) (r > R). (99) 

The electric fields 

E(r) = — Gradf* <p =- ^ ^ e r (100) 

dr 

for potentials (|%ll and (|97l) are 


E(r) = ((£) - 3) C 2 r 2 ~ D + ^ ) e r (0 <r<R), (101) 


E(r) = (D — 3) C 4 r 2 ~ D e r (r > R). 


( 102 ) 


Because the electric field (11001) must be finite at all points, and r 2 ~ D —y oo for r —» 0 and 
2 < D < 3, it is necessary put C 2 = 0. Conveniently potential normalized by the condition 
<£>(0) = 0, then we get C\ = 0. Because there are no surface charges, then the electric field ( 1 1001 ) 
at the surface of the cylinder r = R is continuous, i.e. the derivative of the potential should be 
continuous. The conditions of continuity of the potential and its derivative at r — R give two 
algebraic equations that allow us to determine the remaining two constants C 3 and C 4 by the 
equations 


pR 2 

2 e 0 (r> — l) 


c 3 + c 4 fi 3 - D , 


(103) 


Then we have 

= pR 2 r = PR 0 ' 1 

3 2 £ 0 (D ^ 3) ’ 4 e 0 (D-l)(D- 3)‘ 


(104) 

(105) 
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As a result, the potential is 


<p(t) = 


pr- 


<P(r) = 


2e 0 {D-l) 

pR 2 pR 0 - 1 


(0 < r < R), 


+ 


„ 3 -D 


2e 0 {D-3) e 0 (D-l)(D-3) 

Using (11001) and (11061 - 1 1071) . the electric field has the form 


(r > R). 


(106) 

(107) 


E(r) = 


pr 


£o (D — 1) 


(0 < r < R), 


o R® 1 r 2 ^ 

E (r)= P e r (r > R). 

£o [-is ~ J-J 

For D — 3, we get the well-known results of non-fractal case. 
Equation (j!08j) can be represented in the form 


E(r) = 


pr 


2£o£ e ff,i 


(0 < r < R), 


(108) 

(109) 


( 110 ) 


where £ e ff,in = (D — l)/2 is an effective permittivity of fractal materials. Consider the charge 
per unit length 


n (D-l)/2 R D -1 

t D — P V D-1 — P r( (£) ^ i) /2) ' 

For D = 3, equation fillip gives the value 73 = pn R 2 for non-fractal charge cylinder. 
(Illip equation (1109p can be represented in the form 


( 111 ) 

Using 


E(r) 


_ 1 _ Tp 

2 7T £q £eff,out T 


(r > R), 


where the effective permittivity 


£eff,out 


(p-i) 

2 ^(3-0)72 r((Z? - l)/2) 


( 112 ) 


(113) 


The electric field in the fractal homogeneous charged cylinder is analogous to the non-fractal 
case up to the factor £ e ff,in■ We have a linear dependence on the distance from the cylinder 
axis for 0 < r < R. Electric field outside the fractal charged cylinder differs from non-fractal 
case. For r > R, we have power-law dependence on the distance from the cylinder axis. In 
addition the electric field outside the cylinder is reduced by the effective permittivity £ e ff }OU t- 
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9 Conclusion 


In this paper, differential operators of vector calculus for non-integer dimensional space is 
suggested to describe fractal media and materials in the framework of continuum models. The 
first and second order differential operators for non-integer dimensional space are proposed 
for rotationally covariant scalar and vector functions. We consider some applications for the 
case of spherical and axial symmetries of the fractal material. Elasticity of fractal hollow 
ball and fractal cylindrical pipe, heat distribution in fractal media, and electrostatic field of 
fractal charged cylinder are described to illustrate the suggested approach. In general, we 
can consider not only Erst and second order differential operators for non-integer dimensional 
space. The differential and integral operators of fractional orders can also be considered for 
non-integer dimensional spaces to take into account non-locality of materials. We can note 
that a dimensional continuation of the Riesz fractional integrals and derivatives [8T1 82] to 
generalize differential and integrals of fractional orders for non-integer dimensional space has 
been considered in [S3]. 
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